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Abstract 

In this work we use Monte Carlo simulations to study the phase behavior of spherical caps 
confined between two parallel hard walls separated by a distance H. The particle model consists of 
a hard sphere of diameter a cut off by a plane at a height x, and it is loosely based on mushroom 
cap-shaped particles whose phase behavior was recently studied experimentally [E. K. Riley and 
C. M. Liddell, Langmuir, 26, 11648 (2010)]. The geometry of the particles is characterized by 
the reduced height x* = xl a i such that the model extrapolates between hard spheres for x* ~ > 1 
and infinitely thin hard platelets for x* ~ > 0. Three different particle shapes are investigated: (a) 
three-quarter height spherical caps (x* = 3 / 4 )j ( D ) one- half height spherical caps or hemispheres 
(x* = 1 /2), and (c) one-quarter height spherical caps (x* = V 4 )- These three models are used to 
rationalize the effect of particle shape, obtained by cutting off spheres at different heights, on the 
entropy-driven self-assembly of the particles under strong confinements; i.e., for 1 < H/x < 2.5. As 
H is varied, a sequence of crystal structures are observed, including some having similar symmetry 
as that of the structures observed in confined hard spheres on account of the remaining spherical 
surface in the particles, but with additional features on account of the particle shapes having 
intrinsic anisotropy and orientational degrees of freedom. The x* = 3 A system is found to exhibit 
a phase diagram that is most similar to the one obtained experimentally for the confined mushroom 
cap-shaped colloidal particles under. A qualitative global phase diagram is constructed that helps 
reveal the interrelations among different phases for all the particle shapes and confinements studied. 



I. INTRODUCTION 

Materials that are able to control and modify light propagation are very promising for 
photonic applications.^ Among these, Photonic bandgap (PBG) crystals, which consist of 
materials having different dielectric constants arranged in a periodic structure,^ have been 
the subject of intensive research during the last few decades. PBG crystals can be fabricated 
by exploiting the self-assembly of polymers, block co-polymers, and colloidal particles into 
different morphologies. P In PBG crystals formed by the self-assembly of colloidal particles, 
the geometry of the particles has been envisioned as one of several possible factors that 
can modify, and possibly enhance, the optical properties of the materials.^"^ Clearly, the 
colloidal particles used for this purpose should have specific and well defined shapes and low 
particle-size polydispersity.^ Advances in new routes of synthesis have made possible the 
fabrication of such anisotropic particles in large quantitiesp^HEl so ft h as become possible to 
produce materials with specific properties by tailoring the individual particles. 

Several examples of crystals formed by the self-assembly of anisotropic particles have been 
reported in the literature. Dimers, formed by fusing two lobes of either equal (symmetric) 
or different (asymmetric) diameters, constitute one of the simplest geometries to model 
anisotropic colloids. PBG crystals formed by symmetric ! 10 1 11 1 and asymmetric dimers^^ 
have been reported. The theoretical PGB for this geometry has also been studied, 10 showing 
that the optical properties of the material can change by modifying the symmetry and the 
degree of interpenetration of the lobes. The influence of particle shape on PBG has also 
been studied in a system consisting of spherical silica particles, whose shape can be modified 
into ellipsoids using ion radiation P 2 ^ Assuming that the material of the particles before and 
after the radiation is not modified, the PBG was found to change as a consequence of the 
anisotropic shape and the lattice spacing of the new structure. 

Over the years, geometrical confinement of colloidal particles has been used as a method 
to control their self-assembly into well defined crystals structures.^ For example, in systems 
formed by spherical particles confined between two planar walls or in a wedge, a cascade of 
different phases is observed as a function of the degree of confinement. Experiments^ 2 ** 2 ^ and 
computer simulations^** 2 ^ of spherical particles concur that as the confinement is reduced 
(i.e. as the separation between the confining planes is increased) the crystal structures follow 
the sequence: 1A -)> 13 ->■ 2D ->■ 2ft -> 2A -> 2ft -> . . ., where the symbols A, B, □, ft, 



and CP denote triangular, buckled, square, rhombic, and prismatic lattices, respectively (see 
note in reference 12%]) . Similar research has been carried out by several groups for anisotropic 
particles, under conditions of strong confinement that approached or enacted a behavior 
in (quasi-) two dimensional geometry. Examples of such experimental and computational 
investigations include studies of two-dimensional (monolayer) systems for hard dimers,^^ 
hard sphero cylinders,^ hard polygonsj*2>SZI superdisksp^ and studies of systems confined 
between parallel walls for thin hard-rods^ and octapod-shaped colloids.^ 

Recently, anisotropic particles whose shape resembles the form of a mushroom cap, have 
been the object of several recent studiesM 3 * 19 * 46 " ^ The cross section of an idealized model of 
these particles is show in Figure ma). It has been observed experimentally that under strong 
confinement such mushroom cap-shaped particles exhibit the formation of the following 
sequence of crystal structures as the confinement is reduced^ 1 A r —¥ 1T> — > IS — > 2D — > 
2A r — y 2S, where S denotes a side rotator hexagonal crystal phase, and the subscripts o 
and r refer to the orientations of the particles in each layer; the former describes lattices 
having all the particles oriented in the same direction, while the latter denotes random 
orientations of the particles. Although the crystal lattices formed by mushroom cap-shaped 
and spherical particles might appear similar, important differences can be observed as a 
consequence of the orientation of the particles in the crystals. In the lA r and 2A r crystal 
structures, for example, the particles in each layer are placed in a triangular lattice with 
random orientations that can be either dimple-up or dimple-down, while in the 1S crystal 
phase a bifurcation in the particles orientation (BPO) is observed with the formation of 
alternating stripes of particles oriented dimple-up and dimple-down with either straight or 
zig-zag buckling. The BPO of the alternating stripes in the 123 maximizes the packing 
of the crystal structure] 46 * 49 * As in the case of confined hard spheres, the 1!B transforms 
into a 2D D crystal phase showing also a BPO as the confinement is decreased. In both 
123 „ and 2D the orientation of the particles in each layer is such that the dimples always 
face towards the walls. In the IS and 2S rotator crystal phases, the particles in each layer 
are placed in a triangular lattice with the particles rotating in their lattice positions, but 
keeping the dimples mostly oriented parallel to the walls. As a consequence of the orientation 
and shape of the particles in the IS phase, the hexagonal order is weaker than that in A 
phases. Recently, the same group studied the self-assembly of similar mushroom-cap shaped 
particles synthesized using seed emulsion polymerization. 13 The phase behavior of the new 



particles was found to be slightly different from that previously reported,®! in particular, in 
the buckled monolayer IS* the particles were now mainly oriented parallel to the walls, and 
in the square bilayer 2D* the particles in the bottom layer were oriented parallel to the walls 
while the particles in the upper layer are dimple-up oriented. These results suggest that 
the phases formed by the particles are sensitive to the shape of the particles and, possibly, 
kinetic traps near metastable states. 

Computer simulation can be used to systematically study the effect of particle shape, 
particularly of orientational degrees of freedom, on the self-assembly and phase behavior of 
anisotropic particles. Simulation studies of anisotropic particles having curved surface area 
have mostly focused on bulk behavior. For example, Cinnacchi et al.™ reported the bulk 
phase behavior of contact lens-like particles using Monte Carlo (MC) simulations. For large 
values of the radius of curvature, approaching the limit of infinitely thin hard platelets, the 
system exhibits the formation of a nematic phase at high densities, consistent with the phase 
behavior of hard platelets,^ while for small values of the radius of curvature the system 
shows the formation of clusters with spherical-like shape. The phase behavior of bowl- 
shaped particles (a shape intended to model experimental colloidal particles concurrently 
studied)^ has also been reported based on MC simulations and free energy calculations.^^ 
Simulations of the bowl-shaped model reveal the formation of several crystal and liquid 
crystal structures depending on the thickness of the bowls which ranged from thin hard- 
bowls to hard- hemispheres. 

In this work, we present the phase behavior of an idealized model of mushroom-cap 
shaped colloidal particles, which are modelled as spherical caps interacting only via excluded 
volume, and confined between two parallel hard walls. Using this model we address the effect 
of particle shape (as the fraction of curved surface area varies) on the self-assembly and 
phase behavior for conditions of strong confinement. Specifically, our study encompasses 
three different hard-core particle shapes: (a) three-quarter height spherical caps, (b) one- 
half spherical caps (hemispheres), and (c) one-quarter height spherical caps The rest of the 
paper is organized as follow: In section II the simulation model is presented, in section 
III the results for the three different particles shapes, including some general remarks are 
described. In section IV the main conclusion are summarized, followed by an Appendix 
where a summary of the simulation techniques are described, including the algorithm to 
determine particle-particle and particle-wall overlaps. 



II. MODEL 

The particle model consists of spherical caps (SC), that is, hard spheres of diameter a 
cut-off by a plane at a height \. Particles only experience excluded volume interactions 
and are confined between two parallel hard walls separated by a distance H . A schematic 
representation of the cross section of the particles is depicted in Figure [TVb) . The particle 
model is fully characterized by its reduced height (aspect ratio) x* — x/ a i an d its orientation 
is completely described by the unit vector u along the principal axis of symmetry. Note that 
X* is also equal to the spherical surface fraction defined as the ratio of the curved surface 
area of the particle to the surface area of the complete sphere, while (1 — x*) is the ratio of 
flat to curved surface areas. In order to address the effect of x* on the phase behavior of the 
system, three different values of the aspect ratio x* are analysed: (a) three-quarter height SC 
(x* = 3 / 4 )> (b) one-half height SC or hemispheres (x* — V 2 )> an d (c) one-quarter height SC 
(x* = 1 /4). Throughout this work, structural and thermodynamic properties are reported in 
reduced units, where the pressure, packing fraction, and separation between the walls are 
given by: P* = Pv p /(ksT), rj = Nv p /V, and H* = H/x, where P is the pressure, ks is 
Boltzmann's constant, T is the absolute temperature, N is the total number of particles, 
u r} = 7rx 2 (3cr/2 — x)/3 represents the volume of one particle, and V is the total volume of 
the system. Details of the simulation techniques used to calculate the phase behavior of SC 
are presented in the Appendix. 

III. RESULTS 

A. Spherical caps with aspect ratio x* = 3 / 4 

The phase diagram for SC with x* — 3 / 4 under confinement is shown in Figure pi The 
phase behavior for this system is richer that the one displayed by confined hard-spheres P*"^ 
For the highest possible confinements, corresponding to 1.0 < H* < 1.33 (the upper limit 
corresponds to H = a), the particles translate in a quasi-two dimensional fashion adopting 
only two possible orientations: either face-up or face-down. The first two stable phases ob- 
served at these conditions are a buckled monolayer 1S with BPO and a random orientated 
hexagonal monolayer lA r . Representative configurations are shown in Figures p[a-b). As 
mentioned in the introduction, confined hard spheres also exhibit the formation of buck- 



led and hexagonal monolayers, however, the main difference is that SC can adopt several 
orientations depending of the degree of confinement. 

For 1.35 < H* < 1.53, high-density systems exhibit the formation of a dimer phase ID, in 
which the particles orient parallel to the walls pairing up in dimers that can take on several 
configurations.^ In Figure Km c) one of the dimer crystal phases with a herringbone config- 
uration is shown. Upon expansion the ID crystal melts into a fluid side phase monolayer 
IS in which the particles remain oriented mainly parallel to the walls. For H* = 1.35 (a 
confinement corresponding to hard walls separated by a distance 1.01a"), the system exhibits 
a re-entrant hexagonal ordering behavior upon lowering the concentration: the hexagonal 
order first drops (forming a IS phase), then climbs up and finally drops down again. While 
the IS phase occurs for 1.35 < H* < 1.53, the hexagonal- order re-entrant behavior is only 
apparent for H* values near 1.35. The phase behavior at these conditions is shown in Figure 
[i| where the equation of state (EoS) and hexagonal order parameter \l/ 6 are plotted as a 
function of the packing fraction 77. The only discontinuity observed in the EoS is for the 
ID —7- IS transition. The low-density branch of the EoS, that encompasses the IS phase, de- 
creases monotonically showing no evidence of additional phase transitions. However, Figure 
W[h) reveals that for H* = 1.35 a peak in hexagonal order occurs that spans 0.46 < i] < 0.51. 
For these states, very long runs were carried out to rule out the possibility of artifacts due 
to slow relaxation and mobility; also, additional compression runs started from an isotropic 
state reproduced the same re-entrant behavior. 

To explore the correlation between side orientation and hexagonal order, the fraction 
of particles oriented parallel to the walls f S id e is shown in Figure [41(c). A particle i was 
considered to be oriented parallel to the wall if |u.j • u w | < 0.5, where u w is the unit vector 
perpendicular to the walls; this angle ensures that the rim of the SC curved surface faces the 
walls. It can be observed in Figure EFc) that for 1.35 < H* < 1.53 the fraction of particles 
side-oriented decreases when the system expands. Figure [5] shows snapshots and correspond- 
ing two-dimensional structure factors 5*(q) for H* = 1.35 in the 0.46 < i] < 0.51 range; it is 
observed that when most particles remain oriented parallel to the walls (e.g., at rj = 0.516) 
the resulting anisotropic shape presented for particle-particle interactions frustrates the es- 
tablishment of hexagonal order, while when about half particles orient perpendicularly to 
the walls, they effectively interact as spheres which can enhance the hexagonal order if the 
volume fraction is sufficiently high. It is then a combination of high volume fraction and 



low enough f si d e that promote hexagonal order. 

The peak-behavior (and strong fluctuations) in \l/ 6 values observed in the re-entrant 
region for H* = 1.35 can be explained by its proximity to the lA r phase that occurs 
(in the same 77 range) for H* just below H* = 1.33 (see Figure pi); indeed, re-entrant states 
around H* = 1.35 could be considered to still be "part" of the lA r phase domain. Of course, 
hexagonal order gets accentuated once H* is not large enough (< 1.33) to accommodate side 
orientations hence leading to smaller f si< i e values. Re-entrant behavior notwithstanding, it 
is also difficult to pinpoint the limits of stability of the IS phase, but an approximate phase 
boundary can be sketched based on the f si d e vs. density plot (see Figure |4](c)). A threshold 
of fside = 0.75 is used to draw the approximate phase boundary represented by thick dashed 
lines in Figures [2] and |4} 

When the confinement is decreased, the 12) phase transforms into a square bilayer 2D . 
This ID —7- 2D D transformation seems to occurs with a Xi intermediate, which is formed 
by the separation of some of the dimers in the ID phase, and their reorientation to face 
toward the walls (see Figure |3](e)). As in the case of confined hard spheres, the 2D phase 
transforms into a 2A mediated by a rhombic phase 2"Jl . Likewise 1S , these three phases 
show BPO. As H* increases at intermediate densities, the 2A transforms continuously into 
a hexagonal bilayer with random particle orientations 2A r , which then transforms into a 
hexagonal bilayer 2§ with particles acquiring orientations parallel to the walls. These three 
phases are shown in Figures |3Vg-i). This sequence of transformations 2A —$■ 2A r — > 2S 
seems to occur continuously as no signs of first-order transitions were observed during the 
simulations. 

At high densities, the 2A D phase transforms into a dense dimer side bilayer 2D structure 
(shown in Figure [3tl)) as H* is increased. The 2A — > 2D transformation seems to occur 
through two intermediaries X2 and X3, which are shown in Figures |3Vj-k). In the X2 phase 
the dimers are intercalated between stripes formed by particles with opposite orientations 
facing toward the confining walls. In the second intermediary structure X3, the particles 
forming the dimer in the X2 detach and reorient to adopt a side configuration, while the 
particles facing toward the confining walls (in X2) remain in a similar configuration. Finally, 
the particles that were facing the walls reorient to form the 2D structure. Upon expansion 
at constant confinement H*, the X2 and X3 structures transform into a 2A r . phase, while the 
2D crystal transforms into a 2S phase. These transformation are marked in Figure [6j where 
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the equation of state, bond order parameters ^4 and ^/q, and fraction of particles oriented 
parallel to the walls f S id e are shown for three different confinements: H* = 2.00, H* = 2.40, 
and H* = 2.67. The densest crystal structures for each of these confinements are 2A , 2X 3 , 
and 22), respectively. 

The transformation observed upon expansion are: 2A D -> Fluid for H* = 2.00 (Figure 
ga)), 2X3 -> 2A r -> Fluid for H* = 2.40 (Figure gb)), and 22) -> 2S ->■ Fluid for 
iJ* = 2.67 (Figure gc)). The three structures 2A a , 2A r , and 2S show large values of the 
bond order parameter \l/ 6 revealing their hexagonal order. However, the f S id e profiles are 
different for each of the confinements. For the H* = 2.00 and 77 > 0.56, f S id e vanishes 
indicating that there is a BPO in the structure (see Figure pig)); however, f si d e increases 
once that the fluid region has been reached. For the H* = 2.40 system the f S id e profile 
shows that approximately half of the particles are oriented parallel to the walls, confirming 
that in the hexagonal structures 2A r particles have random orientations as shown in Figure 
gh). Note that while in the lA r phase the particles can take only face-up or face-down 
orientations, in the 2A r they can adopt any possible orientation. Finally, in the H* = 2.67 
system most of the particles are oriented parallel to the walls, as revealed by the large 
values of f S ide, confirming the existence of the 2S phase. It can also be observed that both 
the 2A r and 2S phases show moderate values of the order parameter ^>q, while the 2A 
phase exhibits larger values for this parameter. This effect is mainly caused by the particles 
oriented parallel to the wall perturbing the hexagonal order J*^ 

The SC model with \* — 3 / 4 presented in this section resembles closely the shape of the 
mushroom cap-shaped particles reported in reference IH3 The comparison of the structures 
formed by such SC with those found experimentally reveals that the model is appropriate to 
describe the phase behavior of the real system. As mentioned in the introduction, mushroom 
cap-shaped particles reported in references 13 and 46 exhibit the following crystal sequence: 
lA r -> Yh (or IS*) ->■ IS -> 2D (or 2D*) ->■ 2A r ->■ 2S. Our model of SC with 
X* = 3 / 4 encounters the sequence lA r -> 123 ->• IS ->• 12) — >■ 2D ->■ 2A r ->■ 2S for 
the horizontal band of packing fractions between rj = 0.56 — 0.57 in Figure |2j suggesting 
that the crystal structures in references [13] and H6] are in equilibrium. Figure [7] shows the 
experimental structures for lA r , IS*, IS, 2D , 2A r , and 2S formed by mushroom caps, 
revealing good agreement with those found in simulation (see Fig. [3]) with one exception: in 
the experimental buckled monolayer IS* the particles orient with the dimples parallel to the 



walls, while in the simulated 1S the particles orient perpendicular to the walls. Moreover, 
the simulations reveal six additional phases that have not been observed in experiments, most 
of which occur at higher densities and whose structure may not be accessible to experiments 
due to kinetic trapping. Indeed, the phases observed experimentally appear to be the first 
ordered structures encountered in our phase diagram (for a given H*) while going up in 
concentration. It is remarked, however, that the boundaries of the simulated phases are 
only approximate and (non-trivial) free energy calculation would be needed to refine them 
and to determine the stability of the additional phases. 

B. Spherical caps with x* = V 2 

The phase diagram of confined SC with aspect ratio x* = 1 / 2 i corresponding to hard 
hemispheres, is depicted in Figure |8j As in the previous case, the first stable crystal struc- 
tures observed for this system is a YS> phase, which is now stable for 1.0 < H* < 1.4. Over 
a very small region in the phase diagram, it is also observed the formation of a random 
hexagonal monolayer lA r , characterized by random face-up and face-down orientations of 
the particles. As in the case of SC with x* — 3 / 4 5 this system exhibits the formation of 
three bilayer structures 2D , 2Jl , and 2A. In Figure [9] representative snapshots for these 
phases are presented. For high densities and H* > 1.85, which correspond to wall sepa- 
rations larger than the diameter of the particles, a rectangular plastic bilayer phase 2C P 
is formed, which is characterized by alternating arrays of face-to-face sphere-like dimers, 
and column-like dimers. At high densities the orientations of the sphere-like dimers are 
frozen as shown in Figure [91(e) ; however, when the system is expanded the dimers start to 
freely rotate around their centres of mass as in a plastic rotator phase as seen in Figure 
[9Ff). The 2C p phase resembles the X 2 structure seen for x* — 3 / 4 > i n which the dimers are 
intercalated between stripes formed by particles with opposite orientations facing towards 
the confining wall. However, in the X2 structure the dimers cannot rotate as in the case of 
hemispheres. The 2C P phase also resembles the face-centered cubic rotator phase formed by 
hard-hemispheres at bulk conditions.^"^ 

The phase behavior of hard hemispheres suggest that the stabilization of hexagonal mono- 
layers over large densities and strong confinement is attained only for systems comprised 
of particles with large spherical surface fraction or x* ■ This effect can be explained using 
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Figure [TO] and the following geometrical arguments. If x* is increased, the SC bases (capping 
discs) are pushed apart making the monolayers more stable (see Figure [To|(b)). In this case, 
the distance d between circumferences of nearest discs in alternating stripes increases and 
the contact angle /3 between nearest particles in adjacent stripes reduces: the limit values 
for a perfect monolayer are d = 2a and f3 = 90°. Likewise, it is expected that in SC with 
X < V 2 under strong confinement, the bases of nearest particles in alternating stripes can 
come closer increasing the stability of bilayers over that of monolayers. In particular, it is 
expected that for infinitely thin spherical caps (x* — > 0) the system would form bilayers 
directly without being preceded by monolayers. The same arguments also suggest that the 
stability of the buckled phase should be reduced as x* ~^ 0. These trends are consistent 
with the results presented in the following subsection for x* < V 2 particles. 

C. Spherical caps with x* = V 4 

The phase diagram for SC particles with x* — V 4 under different confinements, shown in 



Figure 11 is simpler than those of the other two particle shapes considered in this work. As 
in the two previous cases, the first crystal phase observed for this system at high densities 
and strong confinement is the buckled monolayer 123 . However, this phase is only observed 
over a very small region of the phase diagram: H* < 1.20 and at high densities. This 
behavior is different from the one observed in confined hard spheres and in the previously 
described SC with x* — V 2 an(1 X* = 3 / 4 > where the buckled phases are stable over a wide 
range of densities and confinements. As a consequence of the small degree of buckling, at 
high densities the 1!B phase shows either straight or zig-zag stripes.^ In further contrast 
with hard spheres, SC particles with x — V 4 do not exhibit the formation of an hexagonal 
monolayer, which confirms the destabilization of monolayers for small x as discussed before 



in reference to Figure 10 The small island of 1!B observed in Figure 11 is expected to get 
reduced if x* is further decreased. For H* < 2.0 the most stable crystal phases observed are 
the 2D and 2 A structures and, as in the case of the previous systems, the phase transition 
2D —¥ 2A is through a rhombic bilayer 2Jl intermediary. All these bilayer structures, as 
in the previous cases, also exhibit a BPO in such a way that the base of the SC always face 
towards the confining walls, allowing them to pack efficiently Representative snapshots for 



the 123 D , 2D D , 2Ji , and 2A phases are shown in Figure 12 
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For H* ~ 2.0 the system shows an abrupt change in phase behavior where the number 
of layers increases from two to four and the structure adopts a square symmetry (4D ). In 
the 4D phase the two middle layers can either align to form dimers as in Figure [l2^e) or 



overlay halfway as in Figure 12 f). In this region, it is difficult to pinpoint which crystal 
structure corresponds to a particular confinement unequivocally. As in the case of the bilayer 
structures, the 4D D if formed by alternating layer of particles with same orientations, with 
the outer layers having only particles oriented facing towards the walls. 

D. Global remarks on the phase behavior for spherical caps 

Our results are summarized in the qualitative x* vs - H* phase diagram presented in 



Figure 13 for the main ordered phases and the 1.0 < H* < 2.4 range only; they are presented 
in the context of the results of the confined hard-sphere system (x* = 1) from reference |2"T1 
Note that hard spheres particles no longer have an orientational degree of freedom and hence 
the phase lA r simply becomes 1A, 123 becomes 12pHl (denoted as 22 and II in reference 
2"7j) . 2CR becomes 2CR, 2D becomes 2D, and phases 2A , 2A r , and 2S fuse into 2A. The 



following trends can be observed: 

• Two phases that occur in all cases are the bilayers 2D and 2A , which are mediated, 
at very high densities only, by the also common phase 2Ji . The lower H* bound in 
the range where the 2D occurs tends to shift to higher values as x* increases, while 
its upper bound remains almost unchanged (at H* ~ 1.8 — 1.9). 

• For H* < 1.5, as x* increases the 2D phase gives way to the triangular monolayer 
lA r (or 1 A for %* = 1). This process seems to always be mediated at high densities by 
the buckled phase 133 (or IS for x* = 1). It appears that the IS and ID monolayers 
observed for x* = 3 / 4 system merge into the 1A and the IS, respectively, as x* 
approaches 1. 
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For 2.0 < H* < 2.4, the 4D D phase seen at x* = V 4 with dimer formation (see Figure 
12^e)) appears to be the precursor of the 2C P phase observed for x* = V 2 , which in turn 
would be the precursor of the X2 phase observed for x* = 3 / 4 5 which (together with X 3 
and at high densities only) would give way to the 27 phase observed for x* = 1. For 
intermediate densities, the 2A r phase observed for x* — 3 / 4 likely originates from the 
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2A phase already occurring for x* = V 2 ! finally, the 2A D and 2A r seen for x* — 3 / 4 
merge into 2 A phase as %* — > 1. 

The 2D and 2 A phases are readily rationalized based on considerations of ordered pack- 
ing of particles that (due to confinement) interact laterally with their rounded peripheries. 
The absence of a lA r . monolayer for the x* — V 4 case occurs because even for H* — > 1 the 
system can accommodate two layers (2D ) by the alternating flipping of particle orientations. 
The 123 D phase can be seen as the best packing of particles at high densities when the system 
confinement is such that neither a monolayer nor a bilayer would fit well. Likewise, the 4D , 
2C P , and X2-X3 phases can be seen as variants of a common motif (where an intermediate 
layer made of "dimers" is sandwiched between two monolayers facing toward the walls) that 
occurs as the system tries to pack densely when the confinement ( 2.0 < H* < 2.4) is such 
that neither bilayers nor higher multilayer arrangements would work well. 

IV. CONCLUSIONS 

The phase behavior for SC under strong confinement inside a slit has been mapped out 
using MC simulations. Depending on the particle aspect ratio and concentration, different 
ordered structures have been observed. The system comprising of SC with x* = V 4 forms 
stable bilayers with square and triangular symmetry, and four-layer structures with square 
symmetry. These systems do not exhibit the formation of stable monolayers, with the 
exception of a very small region where a buckled monolayer is observed. The stabilization 
of the monolayers is only observed when the spherical surface fraction of the particles x* is 
increased. Indeed, stable buckled monolayers are observed over a wide range of densities for 
SC with x* — V 2 (hemispheres) and x* = 3 / 4 > while stable random hexagonal monolayers are 
only observed for the system with x* = 3 / 4 - A common feature of the three particle shapes 
studied is that for confinements where the particles can fully rotate, the system tends to 
form dimers to maximize the packing of the structures. The progression of phases resulting 
from changes in confinement and particle aspect ratio, can be traced and fully integrated 



in a global phase diagram (Figure 13), which further provides a seamless connection to the 
known phase behavior of confined hard spheres. 

The system of SC with x* = 3 / 4 shows a phase behavior similar to the one found ex- 
perimentally for mushroom cap-shaped particles reported in reference 061 including the 
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formation of mono and bilayers with particles having orientations parallel to the walls 
(side phases). The formation of these structures is only possible when the particles 
have a shape close to spherical. The sequence of structures observed for this system is: 
lA r —)■ 123 — > IS — > 2D — > A r — > 2§. Our study reveals that such SC can form other 
structures (at higher densities) that have not been observed in experiments yet. Some of 
the phase boundaries are difficult to pinpoint, and future free energy calculations can help 
elucidate more completely the phase diagram for these systems, and the stability of some of 
the structures. Attempts to obtain the structures by compression runs have met with partial 
success only, revealing the presence of potential kinetic traps. Elucidating the kinetics of 
some of these order-disorder phase transitions using specialized simulation techniques is part 
of our ongoing efforts. 

Appendix 
Simulation details 

The phase diagrams for each particle shape are obtained using MC simulations, and 
consist of two stages. The first stage corresponds to the prediction of the densest crystal 
structures for each pair of parameter values (x*,H*)- The second stage consists of expan- 
sion runs (in the isobaric-isothermal ensemble) from the crystal structures obtained in the 
previous stage, and analysing the different structures observed in the system to assign their 
phase identity. 

The calculation of the crystal structures is carried out using the so-called floppy-box 
Monte Carlo algorithm) 54 * 57 ! In this method, systems comprised of a small number of particles 
are simulated using isothermal-isobaric (NPT) ensemble MC runs, allowing the area and 
shape of the box in the x-y plane to change, keeping the height H* along the z-direction 
constant.^ The methodology entails first the generation of a low density configuration, 
with particles in random positions and orientations. The system is then equilibrated at the 
pressure of P* = 2.5, and subsequently compressed in steps of AP* = 2.5 until the a value 
of P* = 10 is reached, and subsequent values increased by a factor of 10 until a value of 
P* = 10 6 is reached. In some systems, the compression between P* = 10 and 100 was 
carried out using steps of AP* = 10, to allow more time for the particles to reorganize, 
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which is particularly important for confinements close to the boundaries between different 
crystal phases. The need of a gradual compression was also pointed out in reference 1531 for 
the study of the crystal structures of hard-hemispheres at bulk conditions. A set of 50 to 
100 of such compressions with different initial configurations are performed, and the last 
configuration in each compression series is recorded for post-processing analysis. Several 
system sizes were tested in order to ensure that the crystal structures are not affected by 
an insufficient number of particles during the calculation, however it was observed that 
using at least 4 particles per layer was enough to obtain reliable results: larger systems 
(more particles per layer) show similar structures as those with system of 4 particles per 
layer. In the post-processing analysis several order parameters were used to characterize 
the crystal phases. The symmetry of each layer is obtained using bond order parameters, 
defined as: ^ = jjJ2j=i ^~ Sfeii ex P(^ n ^'fc)' f° r n = 4 5 6, and, 8, where 9jk is the angle 
made by the bond between particles j and its nearest neighbor with respect to an arbitrary 
axis, and rij is the total number of nearest neighbors of particle j. The calculation of the 
nearest neighbors is carried out using the Voronoi tesselationp^ except for the case n = 4, 
where only the four closest neighbors are used. The reason of this is to avoid the well- 
known problem of degeneracy in the Voronoi tessalation in square lattices. The square 
symmetry is easily characterized using the order parameters ^s taking into account all the 
nearest neighbors^ however ^ 4 , calculated with the above definition, in combination with 
\l/6 are useful to characterize layers with rhombic (oblique) and rectangular symmetries. In 
addition, the following angular order parameters were calculated: Si = J2i=i cos ($«) an d 
S 2 = X^i=i(| c °s 2 (^i) — §); where 6 1 , is the angle between the orientation of particles ti; 
and the main director. The main director of the system is obtained from the eigenvector 
that corresponds to the largest eigenvalue of the so-called Saupe tensor defined as: Q a ,p = 
^(ui >a Ui^) — |, with a, ft = x, y, zP^ The order parameters Si and S 2 are useful to determine 
whether or not a system shows a bifurcation in the particles orientation. For example S 2 
allows to distinguish between 2A D (large values of S2) and 2A r (low or moderate values 
of 5*2). The structures showing the higher order, measured by the aforementioned order 
parameters, are selected as the crystal structures, which in most of the cases also correspond 
to the densest structures. However, non-trivial free energy calculations may be required to 
unambiguously determine the most stable crystal structures for each confinement, especially 
to pinpoint the boundaries between different phase domains. 
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Once the crystal structures are determined, we proceed to perform expansion runs to 
obtain the phase diagram for each particle shape and confinement. As in the previous 
calculation, the expansion runs were carried out using NPT MC simulations. During the 
simulations involving ordered structures the shape of the box is allowed to change, keeping 
the values of H* constant. The number of particles N used for each system range between 600 
and 700. The system is started at high pressures and subsequently expanded until the system 
is completely isotropic. One MC cycle is defined as N moves consisting of translation moves, 
rotations and inversions of the particle orientations, as well as changes of the volume and/or 
shape of the system box. Each move is randomly selected using the following probabilities: 
35% of translation moves, 35% of rotations, 25% of orientation flips, and 5% of volume 
change (50% of which correspond to isotropic changes and 50% for simultaneous changes 
of the volume and shape of the simulation box). For each state, 2.5 x 10 5 MC cycles are 
used to equilibrate the system, and additional 1 x 10 6 cycles, divided in 10 blocks, are used 
to collect ensemble averages and estimate uncertainties. Once the equilibration run for a 
specific state is finished, that configuration is used as a starting point for the next run at a 
different pressure. The location of phase transitions are identified based on discontinuities 
and inflexion points found in the equation of state, and in the order parameters previously 
defined. 

Overlap algorithm for spherical caps under confinement 

The heart of the MC simulation program is the algorithm to check for the particle- 
particle and particle-wall overlaps. To perform these tests, each particle is modeled as a 
void SC (bowl) capped with a circular disc, thus the particle-particle overlap algorithm 
consists of three elementary tests: bowl-bowl j 53 * 62 ! disc-disc,^ and bowl-disc. The particle- 
wall algorithm can also be separated in two elementary tests: bowl-wall and disc-wall. 



Sketches of the different tests for this algorithm are shown in Figure [14 

The algorithm is started by checking the overlap between the full spherical surfaces of 
particles i and j. If the distance between the centers of the spheres |ry| is larger than the 
diameter of the particle ( | r^- 1 > a) then no overlap between the two particles is possible. If 
|r.y| < a then the elementary tests are considered. 
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1. Bowl-Bowl test 

This test is similar to that in reference [53, and proceeds as follow: 

• Check first if the bowl of particle i overlaps with the sphere of particle j. For this 
purpose calculate the angles Wy and faj, the former corresponding to the angle between 
the orientation of bowl i, represented by the unit vector Uj, and the interparticle vector 
ry, and the latter corresponds to half the opening angle of the cone whose vertex is 
the center of the sphere of particle % and the base is a disc of diameter d c formed by 
the intersection of the spheres of particles i and j. The angle 9 (unique for any particle 
type) is also calculated. These three angles, shown in Figure [T4|( a), are obtained from: 

cos(w ij ) = -^ — p, (A-l) 

I T ij I 

cos(4>ij) = —^-, and (A-2) 

a 

cos(6) = (| - x) \- (A-3) 

The bowl of particle % intersects with the sphere of particle j if 

\cuij + \(f) i:i \ < 6, (A-4) 

for A = 1 or — 1. The intersection between the two spheres of particles i and j is a 



circle of diameter d c . Eq. A-4 holding for both As means that the intersection of such 



a bowl with the sphere of the other particle is also the complete circle of diameter d c 



see Figure |I4[a)). If Eq. A-4 holds only once for one A means that the intersection 



is an arc which is a segment of the circle of diameter d c . If Eq. A-4 does not hold 



for both A = 1 and A = — 1 then no overlap is possible. Otherwise, if Eq. |A-4| holds 



for both As and it holds at least once for the bowl of particle j (performing a similar 



analysis), then an overlap exists. If Eq. A-4 holds once for each bowl, then we proceed 
to the following step. 

Both bowls overlap if the intersecting arcs of each bowl intersect. Otherwise, the 
bowls do not overlap. For this, angles ji, jj, and a^- (see Figure |l4|a)) are calculated 

f rQn |53l62l 
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oos( 7i ) = COS(6,) "!°! ( ^ )CO ? ( ^ ) , and (A-5) 

cos(a^) = , ' I, 3 (A-6) 

I 2 11 J I 

where u+ is a vector in the plane of the intersection disc that cuts in half the intersec- 
tion arc of particle i, and is obtained as: 

uf = Ui- {Tij ' f )Tij . (A-7) 

The arcs of bowl % and j overlap if 

M < N + IT? I- (A-8) 

2. Disc-Disc test 

For a spherical cap of diameter a and aspect ratio x> the radius of the capping disc r^ 
is given by: r^ = (o" 2 /4 — h 2 ) 1 ^ 2 (see Figure 14 b)). The distance between the center of 



the discs of particles i and j is r^ = r^ + (a/2 — x)(iij — Uj). Following reference [63], the 
disc-disc overlap test is carried out as follow: 

• If the distance between the centers of discs i and j is larger than the diameter of the 



discs, |r^-| > 2r^, then no overlap is possible (see Figure 14 'h)), otherwise the test is 
continued. 

• Calculate if the plane belonging to the disc i intersects with the plane of j and vice 
versa. If either intersection does not exist, then no overlap is possible. This test if 
performed calculating the distance a, from disc i to the intersection of the planes of 
discs i and j, and it is obtained as 

(r d ■ u-) 2 
^ = A ', L , (A-9) 



and likewise for disc j: 



l-(u,-u^ 
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' _ l-(Vfi)) 2 ' ( °' 

Disc i intersects the plane of disc j if a 2 < r; 2 ,, and so does disc j if a 2 < r 2 ,. 

If both discs intersect the planes of each other then we proceed to check if the discs 
overlap. Overlaps between the two discs exists if the following inequality holds: 




, dj ^\6 tj ) + (r i y-a i )(u i -u j )(r i ^u j ) (A-ll) 

< [(r 2 sin 2 (^) - (r^ • u,) 2 ) (r 2 sin 2 (<y - (r^ ■ u 3 ) 2 )] 1/2 
where 8^ is the angle between the vectors iij and iij. 

3. Bowl-Disc test 

The bowl-disc overlap test is carried as follow: 

• The shortest distance d from the center of the sphere of bowl % to the plane containing 
the capping disc of particle j is calculated (see Figure [l4Fc)). This distance is given 



by: d = \r t b - ■ Uj|, where r t b - = r^ + (er/2 — x)uj is the vector from the center of the 
sphere i to the center of disc j. If d > a/2 then no overlap is possible. If the d < a/2 
then we proceed to the following step. 

• The intersection between the plane containing the disc of particle j and the sphere of 
particle i is a circle of radius r a = ((a/2) 2 — d 2 ) 1 ^ 2 . We call p the shortest distance 
from the center of disc j to the center of the intersection circle of radius r a , which is 
denoted as d and is calculated as p = (|r.;j| 2 — d 2 ) 2 . If (p — r a ) > r d then no overlap 
is possible. 

Note that this test can be optimized by excluding some orientations that also include the 
bowl-bowl and disc-disc overlap tests. More specifically, if 

cos(0) < cos(<£y) = iii ■ iij < 1 (A-12) 

then the bowl-disc overlap test can be avoided as this implies that disc-disc and/or bowl-bowl 
can occur along the bowl-disc overlap. 
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4- Bowl-Wall and Disc-Wall tests 

Both the bowl-wall and disc-wall overlap can be considered simultaneously in the following 

way: 

• We computer the angle /3j between the orientation of particle % and the unit normal 



vector perpendicular to the wall u w (see Figure 14 d). This angle is obtained as 
cos(/3j) = Uj • u w . If cos(/3) < cos(#) then the only check that has to be performed is 
the one involving the sphere of particle i with the plane of the wall, i.e., if the distance 
d w from the center of the sphere to the wall is d w > a/2 then no overlap is possible. 

• If cos(/3j) > cos(#) and the center of the capping disc lays outside the confinement, 
then an overlap exists. 

• Check if the disc of particle i and the plane of the wall overlap in analogy with with 
the disc-disc overlap test. 
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(a) 



(b) 




x* = 3 A x* = l fr x* = V4 

FIG. 1: (a) Cross section of an idealized model of mushroom cap-shaped particle taken from 
reference H6l (b) Cross section and geometry of the spherical caps particle model studied in this 
work. The particle model consists of a sphere of diameter a cut off by a plane at a height x- 
The geometry of the particles is characterized by the reduced height (aspect ratio) x* = x/ a -> an d 
the orientation is defined by the unit vector u along the principal symmetry axis of the particle. 
Depending on the aspect ratio, three different models are defined: (c) three-quarter height spherical 
caps (x* = 3 / 4 )> (d) one-half spherical caps or hemispheres (x* = 1 /2), and (e) one-quarter height 
spherical caps {x* = 1 /4). For each model the bottom, side and top views are shown. 
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FIG. 2: rj-H* phase diagram for spherical caps with x* = 3 / 4 obtained from expansion runs using 
NPT MC simulations. Dashed lines represent approximate phase boundaries. 
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FIG. 3: Representative snapshots for the different phases formed by spherical caps with \* = 3 / 4 
under confinement, (a) Buckled monolayer 123 , (b) random hexagonal monolayer lA r , (c) dimer 
monolayer 12), (d) side phase monolayer IS, (e) intermediate Xi monolayer, (f) square bilayer 2D , 
(g) hexagonal bilayer 2A D , (h) random hexagonal bilayer 2A r , (i) hexagonal side phase bilayer 
2S, (j) and (k) correspond to the intermediate X2 and X3 bilayers, and (1) dimer bilayer 2D. The 
particles are coloured according to their orientation with respect to the axis perpendicular to the 
walls. Only a small section of each configuration is shown for clarity. 
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FIG. 4: (a) Equation of state in the P* — r/ plane for spherical caps with \* = 3 / 4 under different 
confinements: H* =1.35, 1.40, 1.47, and 1.53. (b) Hexagonal order parameter \&6 as a function 
of the packing fraction r/ for the same confinements shown in (a), (c) Profiles of the fraction of 
spherical caps oriented parallel to the walls f s ide as a function of r\ for the same confinements 
shown in (a) and (b). The lines are polynomial regression fits of the simulation points. The 
dotted (brown) lines shown in both (a) and (c) denote the approximate boundaries of the IS phase 
corresponding to states having 75% of their particles oriented parallel to the wall, i.e., points laying 
on the constant line of f s ide = 0.75 in (c). The points in (a), (b) and (c) marked with a large (red) 
open circle correspond to the states shown in Figure [5J 
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FIG. 5: Configuration snapshots and two-dimensional structure factors S*(q) for spherical caps 
with x* = 3 / 4 an d H* = 1.35. The states correspond to the points marked with large (red) 
open circles in Figure Hi From left to right these states are: P* = 12.1(r? = 0.516, l^ = 
0.105, f side = 0.819), P* = 11.1(7? = 0.505,1-6 = 0.411, f side = 0.753), P* = 8.0(77 = 0.473, ^ 6 = 
0.702, f side = 0.517), and P* = 6.6(7? = 0.454, $ 6 = 0.165, f side = 0.517). S(q) is calculated as 
S(<i) = (1/^0 (E«=i cos(q • Yi)] 2 + Ei=i s i n (q • rj)] 2 ), where rj is the position of the centroid of 
particle i and q is a reciprocal wave vector. 
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FIG. 6: Equation of state for spherical caps with \* = 3 / 4 f° r three different confinements: (a) 
H* = 2.0, (b) H* = 2.40, and (c) H* = 2.67. For each confinement the pressure P* , order 
parameters ^4 and *$?$, and fraction of particles oriented parallel to the walls f s ide are shown as a 
function of the packing fraction 77. 
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FIG. 7: Experimental mesophases based on the self-assembly of polystyrene "mushroom caps" 
under wedge-shaped confinement, (a) lA r , (b) IS*, (c) IS, (d) 2D , (e) 2A r , and (f) 2S. The de- 
scription of the particle shapes and confocal microscopy videos of assemblies are found in references 
[Ljandliej 
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FIG. 8: rj-H* phase diagram for spherical caps with x* = V 2 (hemispheres) obtained from 
expansion runs using NPT MC simulations. Dashed lines represent approximate phase boundaries. 
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FIG. 9: Representative snapshots for the different phases formed by spherical caps with x* = V 2 
(hemispheres) under different confinements, (a) Buckled monolayer IS,,, (b) random hexagonal 
monolayer lA r , (c) square bilayer 2D , (d) hexagonal bilayer 2A D , and rectangular bilayers 2C P 
with dimers in (e) frozen state and (f) plastic rotator state. The particles are coloured according 
to their orientation with respect to the axis perpendicular to the walls. Only a small section of 
each configuration is shown for clarity. 
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FIG. 10: Buckled crystal structures for maximum confinement of spherical caps with different 
aspect ratios: (a) x* = V 4 ) ( D ) X* = V 2 ) ( c ) X* = 3 / 4 > an d (d) X* = V 95 - & denotes the contact 
angle between nearest particles in adjacent stripes and d is the distance between the circumferences 
of nearest discs in alternating stripes. The formation of an hexagonal monolayer is obtained when 
d — > 2cr and (3 — > 90° as x* —> 1, which corresponds to the hard-sphere limit. Conversely, the 
formation of a square bilayer is obtained when d — > and /3 — >• 180° as x* ~^ 0) which corresponds 
to the limit of infinitely thin spherical caps. 
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FIG. 11: rj-H* phase diagram for spherical caps with x* = V 4 obtained from expansion runs using 
NPT MC simulations. Dashed lines represent approximate phase boundaries. 
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FIG. 12: Representative snapshots for the different phases formed by spherical caps with x* = V 4 
under different confinements, (a) Buckled monolayer 1S D , (b) square bilayer 2D , (c) rhombic 
bilayer 23i , (d) hexagonal bilayer 2A , and four-layer square phase 4D with (e) dimers formation 
and (f) overlaying alignment. The particles are coloured according to their orientation with respect 
to the axis perpendicular to the walls. Only a small section of each configuration is shown for clarity. 
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FIG. 13: Qualitative summary phase diagram of main ordered phases encountered in spherical 
caps with different x* values (x* = 1 corresponds to hard spheres) and different confinements 
(H*). For any given x* i each colored bar represents a phase with the bar appearing lower (higher) 
corresponding to the ordered phase occurring at lower (higher) volume fractions. The data for 
X* = 1 is taken from reference 1271 
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FIG. 14: Description of the overlap algorithm for spherical segments under confinement. To 
implement this algorithm, each particle is modeled as void spherical cap (bowl) capped with a 
circular disc, hence the algorithm comprises the following elementary overlap tests: (a) bowl- 
bowlpMli q^ disc-discj^ (c) bowl-disc, and (d) bowl-wall and disc-wall. 
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